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ABSTRACT 


This paper is an application of the Piecewise-Sinusoidal 
Reaction Matching Technique (PSRMT) to linear, parallel 
thin-wire dipoles. The "method of moments" if briefly 
Suciined amd the fundamental electromagnetic principles 
Direc iprocicy, Surface equivalence and reaction are dis- 
@eesced, Phrough a logical combination of these concepts, 
mies developed and applied to thin-wire dipoles radiating 
in free space due to a delta function voltage generator at 
Maeir center. The driving point impedance and the current 
Mrouribpuvion are calculated and compared to results obtained 
from other independent theoretical techniques. 

iiemilaTor Colliers pbulTtonm Of this anvestigation is not that 
Ore presenting a new technique for oaLouiletitne the input im- 
meaance and current distributions on dipole antennas. The 
Seudy 15 primarily concerned with comparing other well-known 
Pechniques with that of PSRMT. Emphasis is placed on the 
ease and convenience of expressing the mutual impedance 


macerrals of coplanar dipoles in closed form. 
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I. INTRODUCTION 


A. DESCRIPTION OF ANTENNA PROBLEM UNDER INVESTIGATION 

tite, paper 15 an investigation of the efficiency of 
the Piecewise-Sinusoidal Reaction Matching Technique (PSRMT). 
In the solution of problems using the "method of moments" 
[3], computation time on the computer depends primarily 
upon the amount of computational effort required to generate 
the elements in the generalized impedance matrix. Usually 
Chis relates to the amount of numerical integration required 
morcalculate the Zon = However, by employing piecewise- 
Sinusoidal basis functions to obtain the impedance matrix, 
Siwesed form expressions can be obtained for the integrations. 
Also, the size of the impedance matrix required to adequately 
solve the antenna problem is relatively small when compared 
GO the use of other basis and weighting functions. For 
wivese reasons PSRMT will be applied and the results compared 
to other techniques. 

The antenna configuration under investigation will be 
restricted to a linear, thin-wire dipole. The antenna is 
werwned tombe radiating invfiree space and is driven at the 
meewer DY a delta function voltage generator. In this 
Peececularvapplication of PSRMT, all segments of the fila- 
mentary dipoles are parallel and of equal length. This is 


Meummemnecessary condition; but one of convenience. The 





Sevieepeine inpedance and the current distributions 


VeemeeDe Calculated DY an application of PSRMT. 


B. BRIEF HISTORY OF RELATED WORK 

In the past much work has been done in the area of 
determining the mutual impedance between two elements with 
paeaceumed Sinusoidal current distribution. Carter [5] 
has presented the mutual impedance between two wires in 
echelon where each antenna is an odd number of half-wave- 
lengths long. The mutual impedance expression between 
two identical nonstaggered parallel center-fed antennas 
of Brown [1] is of historical renown. The work of Brown 
was carried one step further by Cox [6] who developed the 
mutual impedance between parallel antennas of unequal lengths. 

In 1957, H. E. King published a paper [10] that developed 
an expression for the mutual impedance between two staggered 
meralilel, center-fed, infinitely thin antennas of unequal 
lengths. Lewin [12] analyzed the mutual impedance of two 
coplanar half-wave dipoles, in arbitrary positions. Finally 
Lewin's work was generalized by Richmond and Geary [10] 
mom any length antenna. 

MuEEOumeMC above Suudies were restricted to assumed 
Sinusoidal dipole current distributions. Without assuming 
a current distribution upon an antenna the current distri- 
bution may be found by breaking the antenna up into several 
segments, calculating the mutual impedance between each 


pair of elements in turn and generating an impedance matrix 


10 





momeeeeciiaracverlZcs the current distribution. This calcula- 
Gion of the mutual impedance is done in general by the 
exepresston published by Richmond and Geary. The technique 


erucely described is that of reaction matching. 


fe SCOPE AND LIMITATIONS OF THE STUDY 

This paper develops the multisegment reaction technique 
meen fundamental electromagnetic theory. Although the 
technique is applied only to linear thin-wire dipoles of 
eoroce COnNduCctiVity, the procedure may be extended to 
fore general and intricate geometries. To cite one of 
the many examples, PSRMT has been successfully employed 
by Travieso-Diaz [20] to the analysis of antennas and 


meavtverers with arbitrary surface shape. 


ugh 





Pee Owl OMe lie METHOD OF MOMENTS 


ltemuolbowine Giscussi@m portrays the major steps 
puplLoyea in the method of moments. The object is to 
meauce the functional equations of field GE OIela Me oils e oe 
equations which can be readily solved. The name "method 
of moments" has been given to the mathematical procedure 


for obtaining the matrix equations [8]. 


A. FORMULATION OF THE OPERATOR EQUATION 

Consider only equations of the inhomogenous type 
L(f) = g, where L is a given operator, g is the excitation 
or source (known function), and f is the field or response 
(unknown to be determined). 

One must determine from physical considerations that 
meee oroolem an question is deterministic, that is, if the 
solution for f when L and g are given is unique. The method 
of moments gives a peneral procedure for treating field 
problems, but Memee ead ts or solutdion vary widely with 
the particular pecopolemesoihe analysis to follow is”that 


which is germane to the antenna problem under investigation. 


Be POSTULATION OF THE BASIS FUNCTION 


Consider the deterministic equation 


MGs) = 5. (1) 


I 





where L is a linear operator, gis known, and 7s Pompe 
determined. Let f be expanded in a series of functions or 


a linear combination of a set of N vectors {Ft 
mac. fF , (2) 


were the c's are constants that are not known but are 

to be determined in the solution process. The is are 
chosen to be linearly independent and are called expansion 
evriGOas Or basis functions. Wor exact solutions, Equation 
(2) is usually an infinite summation and js form a complete 
meweotr basis functions. For approximate solutions, Equation 
2) is usually a finite summation. Sroeeteucing Equation 

(2) into Equation (1), and using the linearity property 


en i, 
Ze 2 (3) 


C. POSTULATION OF THE WEIGHTING FUNCTION 
It is assumed at this point that a suitable inner 
product <“f,g> has been determined for the problem. The 


moner product 1s a scalar defined to satisfy 


ieee a ¢g,f7 
CK + eg,h> = X<f,h> + 6<g,h> 


Sooo if f 7/0 
[HO f= 0 , 


6 





where & and @ are scalars, and * denotes complex conjugate 
fol The appropriate inner product for this particular 
Mmieoina proolem will be developed in the section of this 
paper entitled "Reaction Concept". Now define a set of M 
ipgearliy independent weighting or testing vectors iw}, and 
take the inner product of Equation (3) with each W » te 


result is 
f. ec, Ww sLf > = <w_,e> (4) 


fe), 2, 3, --- M, where the N c's are unknown and all 


other terms are known [3]. 


D. SOLUTION OF THE LINEAR ALGEBRAIC SYSTEM 


Equation (4) can be written in matrix form as 


lll ® 


where 


mn 


14 





| | C5 a, Wo 98 
C = and g == 
n n ae 


PcouorSvenumequabilon results by requiring M and N to 
Be equal. That is, the basis set ic and the testing set 
We each has N elements. The coefficient matrix is deter- 


mined by inversion: 


HF ‘ 


Mmeevided that the matrix [Zan is nonsingular. The 


Solution for f is given by Equation (2): 


r= B fe, | (7) 


This solution may be exact or approximate, depending upon 
the choice of i and W and the extent of the summation. 


The particular choice We = oe is known as Galerkin's method. 


E. COMMENTARY ON THE METHOD OF MOMENTS 
One of the most critical tasks in any particular problem 
Hsecvne choice of the basis and testing functions, f. and 


1g) 


We The fi's should be linearly independent and chosen so 


IS 





that some superposition (see Equation (2)) can approximate 
f reasonably well. The ws should also be linearly 
independent and chosen so that the products <w 87 

depend on relatively independent properties of g. Obviously 
additional considerations affecting the choice of ie and 

W, are a.) the accuracy of solution desired, b.) the ease 
of evaluation of the matrix elements, c.) the size of the 


matrix that can be inverted, and ad.) the realization of 


a well—conditioned matrix |E,.| C8). 


FF, PIECEWISE-SINUSOIDAL BASIS SETS 

The physical problem may be represented by many 
menereny. operator equations, and a suitable one must be 
@aoaen. There are also an infinite number of sets of 
expansion functions ce and testing functions Wa that may 
wemenosen. 

Miviueivelvowerne more Nearly the expansion set resembles 
the sevual Torm of the eurrent distribution, the fewer the 
terms required to-accurately represent the current. (In 
addition, certain types of basis functions may lessen many 
time-consuming, numerical operations incorporated in 
approximations. The piecewise-sinusoidal basis set exhibits 
these qualities. The objective is to represent the current 
distribution, over individual subsections of a domain with 
€xpansion functions that possess greater curve-fitting 
potential than do flat pulses and yet lead to well-conditioned 


matrices in approximate solution procedures. The application 


16 
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Of piecewise-sinusoids in Galerkin's method and the 
method of moments to electromagnetic field problems is 


known as the reaction concept [3,17]. 


ie 








iti. THEORETICAL MOTIVATION 


The following discussion will develop the basie concept 
ieee Zzed 1n the formulation of the reaction matching 
meennrgque Of J. H. Richmond. The reaction concept of 
Rumsey [17], the reciprocity principle [9] and the surface 
equivalence theorem [18] were presented by Richards [13] and 
Travieso-Diaz [20]. The concepts of reciprocity, reaction 
and surface equivalence form the theoretical motivation 
more PoRMT; therefore they are included again here. With 
feet ical combination of these elements, one can develop 
a general formulation of electromagnetic boundary-value 


maemplems, in particular, reaction matching. 


fee neCLPROCITY THEOREMS 

Consider a pair of time-harmonic sources of the same 
meequency WM radiating in a linear, isotropic, but not 
meeessarily homogenous medium with permittivity ¢€ and 
permeability p. One of the sources, designated as source 
ieeias impressed electric and magnetic current densities 
(J,, M,) which radiate fields (E,, H,) in the presence of 
peeoecond source, source 2, acting as a receiver. Similarly, 
the current densities (F,, ) of source 2 set up fields 
(E,; H, ) in the presence of source 1 acting as a receiver. 

The fields of these sources at every point in the 
medium are related by the point reciprocity theorem of 


Merentz |9| as follows: 
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ii 
——* 
Cy 
[Ti] 

j 
Cy 
[| 
Nee? 


H, ) 
(8) 


| 
=o" 


lemmueuecMecmoources be Of LTinite extent, and assume all 
inhomogeneities in the medium to have a finite distribution. 


Memane non-overlapping surfaces S. and S5 enclosing the 


al 
source distributions (J,; M,) and (Jo, My) respectively; 


Eee inhnomogeneities are contained in the volumes Vy and V5 


1 and So. Mherefore, the medium outside 5s 


and S5 fPmeomeorctcouc me unMacr these conditions, the 


bounded by S 1 
reciprocity theorem of Carson [4] establishes a relation 
between the source distributions (J). M, ) and (J; M,) 


and the fields (E,, H) and (E,, 


Hy) of Equation (8) as 


SS (J, + EB, - Mi, + B)av = Me. /E, - M+ Hav 
1 | : (9) 


B. SURFACE EQUIVALENCE PRINCIPLE 

The surface equivalence principle can be briefly 
Gescribed as follows [13]. First, an original source 
distribution (J, , M, ) set up fields (E, H) everywhere in 
some defined medium. Next, the original source is turned 
off, but simultaneously surface current densities (J. Mm.) 


are set up on a surface S which encloses the original 


IL 





: 
| 





source and any inhomogeneities. If the surface current 
densities (J. Mo) are related to the original fields 


eae H) by 

pan A nee 

J =nxXd (10A) 
== -— “A 

M=EXn_, (10B) 


A 
Weere nm is the unit normal at 5S directed outward from S, 
Bien cChese suriace current distributions will radiate the 
original fields (E, H) in the region exterior to the surface 


© and will generate a null field inside 5S. 


eee RaACTLION CONCEPT 
In 1954 Vv. H. Rumsey published a paper [17] in which 
he introduced a new physical observable, "reaction". 


Define a complex number denoted by <a, b> as follows: 


~~ 


<a, b> = Wftze - E(b) - M(a) + H(b) ]dv, 
V 
| Ca) 


where the volume V contains the source a. 


The reciprocity theorem states that 
ga, b> = <b, ad (12) 


provided that all media are isotropic and that a and Db can 


be contained in a finite volume. The scalar <a, b> is 
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a measure of the reaction (or coupling) between the sources 
Peenioebee fouation (11) is to be understood as a formula 

for the measure of reaction. It is obvious that Equation (9) 
may be restated as Equation (12), where the inner product 
represents the complex integrals of Equation (9). Equation 
(12) is, therefore, the reaction form of reciprocity. 

These theorems will now be used to develop the "Piecewise 
Sanusoidal Reaction Matching Technique" (PSRMT) of J. H. 


Rachmond. 


ail 





IV. PIECEWISE-SINUSOIDAL REACTION MATCHING TECHNIQUE (PSRMT) 


fee LNTRODUCTION TO PSRMT 

In this section, the method of moments and the reaction 
eoncept will be applied to the solution of the current 
Oistribution on linear antennas. This particular class 
eeeeoroblem involves only antennas which consist of or can 
be represented by a collection of thin-wire segments. 

An axial filamentary test source will be developed. 
Aiter specifying a suitable current distribution on both 
Mrpemorifina: ana the vest antenna, the reaction concept 
will be applied. Reaction will be employed in order to 
Gevelop a system of linear- equations whose sclution will 
yield the complex coefficients that complete the description 
meee approximate current distribution specified. 

mimes react1on Concept and its applications have been 
discussed by Rumsey [17] and Harrington [9]. The interpre- 
magion of PSRMT to follow is primarily that of Richards [13] 
ana Leonard Tsai [21]. Most of the material presented 
here has been previously presented by Richmond [14, 15, 16] 


and Richards; it is presented here for completeness. 


B. THIN-WIRE PROBLEM 


Consider the thin-wire dipole shown in Figure 1. The 


1 
necessarily equal, but shorter than 4/2. The radius 


Baole arms are Wire segments of length L. and Ly» not 


| 


a 


ae 








e— Magnetic Current Ring 


G 


hyviel Test Source 


<— ao 


n-wire dipole with @ filamentary 


FIGUrE 1. A thi 
test sourcé cn its @7is. 


nr 
Lad 








of the wire is arbitrary, but in the thin-wire case always 


eieieetin terms of A; also, a << L. and a << L,. The 


1 
feeole 18 driven at the junction of its arms with a magnetic 


mae Current Ma 


bet the dipole depicted in Figure 1 radiate in free 


space, setting up time-harmonic field intensities (Eas Ha) 
everywhere. By the surface equivalence principle this 
maeenna may be replaced by a set of surface currents 
(Jas M5) enclosing the primary source where 
TJ, = nx (13A) 
oil a d 3 
— —_ Aw 
M4 = Eg a ih (13B) 


peeeiene surface, and n is the unit normal vector pointing 
@euward from the surface. Let this surface exist at an 
infinitesimal distance outwardly away from the actual 
emoole, SO that the dipole itself is replaced by surface 
currents (ees Then it is known that exterior to this 


errace the two fields (E Ha) exist and there is a null 


aie. 
field eeice Mee suri ace. Now tne Original Source can be 
removed without disturbing the fields, so that there is 
free space everywhere. | 

Peebest source tS placed at the axis of the actual 


dipole in Figure 1, in order to employ the reaction test. 


Alone, this filamentary test dipole radiates the fields 
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(B,, Fy) 


Equation (9) may be applied between the test dipole and 


MiMi ces alco wine Carson reciprocity theorem of 


Mae Unknown equivalent surface current on the original 


seepole to give: 


Se. weet Helds = Jc, 1S Salih 03 Thee 
: (14) 


where - Towne wsuriace OL Lhe Original dipole as shown in 
Figure 1, £ is the coordinate along the filamentary test 
mipole, and (T,, M,) represents the currents of the test 
meoole. Due to the null field generated by the original 
memrece in the interior of Se: where the coordinate along 
mies tilamentary test dipole is, the right-hand side of 
Equation (14) is identically zero. Therefore the reaction 


Meech Lhe two sources iS Zero, 
a, be = Cb, arm =O , 


hence the term "Zero reaction condition". This zero 
reaction theorem was developed by Rumsey [17]. Equation (14) 


feeeeoen in light of the zero reaction condition is: 


— 
Qu. 
| 
a 
Q, 
62) 
ev) 
] 
(A 
— 
= 
QQ. 
to] 
iz 
Q, 
2) 
ay 


: —_— (15) 


Hor a typical antenna problem, the excitation consists 


of a voltage generator driving the dipole at the junction 


25 





oo ulS arms. 


voltage generator is assumed. 


densities are: 


S 
tl 


where Va 


is the radius, 2 


p Ci ay ek) ; 


Pome Men otic Olumol Une VOlUape Menerator , 


Since the antenna is thin, a "delta gap" 


Here the impressed current 


(16A) 


(16B) 


itech! 
a 


0 is the value of the coordinate 2 at the 


pap and § is the Dirac delta function. 


In this case M 


terminals, 


Equation (15) reduces to a line integral. 


A Pomoc @ Maenel le Cumrent at the 


and the integral on the left-hand side of 


Applying Ampere's 


Seow tO perform the integration around the gap: 


Combining Equations (15) and 


where Wy (£0) is the complex 


dipole. 


rae 
= =a. Ve al H, (a, £5 f)a# 
Cry) 
Gi ay ead. 


SESE ) (18) 


sree ag 


0 


terminal current on the test 
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For the evaluation of Equation (18) allow the thin-wire 
Meproxkimation (a << *) to hold, thus the distribution of 
current density may be assumed to be uniform around the 
merimnder. ihis assumption may be expressed mathematically 
as simply, J = I,/2na where "a" is the radius and Ta 
foeees Only along the length of the dipole. Thus Equation 


(18) becomes 


wae (L,) Ey dl. , (19) 


c 


ey 
tH 
ct |e 


Equation (19) is the result of replacing the surface 


emrrent density q with an equivalent filamentary current 


of strength I, = etrad 5 


Pet ace where Jog Residence Une anlegration dk. is along 


which 1s located on the previous 


this line on the dipole surface and Es represents the field 
Seeeche test dipole on this line. It should be noted that 
meme CFLest dipole were not axially located, then this 
meelification would not be judicious. 

Equation (19) forms view oOlmdayrTom Upon whieh the two 
segment reaction matching technique is based. Once Va 
is specified and Chee Meron eais tribution Tg ew ssc Cs 
the coefficients for the current i may be readily Cetermined. 
For the remainder of eee BoewciscuSsiomewi ll “involve 
Cement distributions on both the test dipole and on the 


original dipole that are "piecewise-sinusoidal"” in form. 


This will lead to considerable simplification because 


e7 





a.) closed form expressions for the Be 


and b.) closed form expressions for [te ° Ey dt . also 


exist for coplanar dipoles. The basis for these claims 


VemmMucaa be found . 


will be presented after a generalization is made for N 
eeenent reaction matching. | 

An essential feature of any acceptable test source 
(J. M,) must be that the fields it radiates be known and, 
df possible, that they be known in closed form. One such 
Beurce 15s a filamentary electrical source with time-harmonic 
current, as shown in Figure e. 

bet tne source be located on the 2 axis with endpoints 


Ze and Zo > Cumvlomevirenued seribubion be of Che form: 


1 


TEC) ON ec) rt (20) 


where A and B are complex constants, k = 2m/a, and the 


ouwt 


time dependence Homer SvOod aad SUDDressed, 


C. BLECTRIC FIELDS OF A PIECEWISE-SINUSOIDAL FILAMENTARY 
SOURCE 
As expression for the Z component of the free-space 
eectoric field of the filamentary source of Figure e is 


given by Schelkunoff and Friis [19] as 


Ho i T! el — 4h et tho a QV (20 
7 aj w il Ry 2 Ry Zo 
— 'es ~ = 
bere I, = I(z,), eae pal = Lz), 1s nC rer 





qe ee _——— 








Figure 2. 


origin 


Filamentary electric line source on the Zz axis 
sets up a field ec! at an observation 
point (@,2). 
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ead V is the electric potential of the end charges; €& 
meetcne permittivity of the medium. Throughout this paper 
mice primne denotes taking the derivative with respect To Zz. 

Mie electric field i is composed of both a normal 
component Ee and a tangential component Eo as depicted 
Meer resure ec. 

Note that in performing the dot product Nd Ei," 
of Equation (18) for linear antennas, only the z component 
@pechne Clectric field makes a contribution. That is the 
ieee rtication for ignoring the Ee component when considering 
eniy linear segment configurations. 

The term 9V/3z will be dropped immediately since it 
vanishes for adjacent segments and the contributions from 
mes pradient of Che potential V due to end charges for the 
end segments may be assumed negligible [13]. Note that 
mmemouch this expression for ES inne abave Ly Sampie, this 
closed form expression is rigorous even in the near-zone. 

The complex coefficients A and B in the expression 
for I(z) (see Equation 20) can be solved in terms of the 
end currents I, and I,. The resulting equation for I(z) 


ul o 
is 


jm) = [it 


, Sin k(Z. ~ 7) + I, sin k(z - z,)J Visio <6! 


(225 


Woere d as indicated in Figure 3 is the length of the 


Seement. 
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FIGURE 3. Piecewise-sinusoidal current distribution 
on an antenna segment. 
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Figure 3 illustrates the current distribution expressed 
mmeomwaction (22). it is clear that from Equation (22) 


I! and I} may also be expressed in terms of I. and I... 


Ik 2 1 2 
Since 
m7) = flat) cos kK(Z5 ~ Zz) + kI, cos k(z - z4)3J /sinkd 
(23) 
imeet1OlLlows that 
t = os 
Ii k (I, I, goon 7 sin kd CA aN 
Is = k(T, eos kd + I,) Pac inekaoy - (24B) 


Note that from Equations (21) and (24) the electric field 
Mimemsci ty Ol the sinusoidal electric Line source may be 
expressed entirely in terms of the geometry of the observation 


Deant and the end currents I, and I,. 


D. EXPANSION OF THE CURRENT DISTRIBUTION 

The forms of the current distributions on both of the 
filamentary Be noiles is piecewise-sinusoidal. The current 
on each segment of each dipole can be described by the 
sinusoidal functions given in Equation (22). Because the 
meeonents of the electric field intensity of such an 
element are known, both terms within the integral of Equation 


(19) exist as closed form expressions. 


Sie 





Mhe following piecewise-sinusoidal distribution 


has been employed by Richmond [15]: 


P. N-1 sin k(2-2£ ) 
5 se ee ee aF5 os Fg Ve n= 
n=1 sin kd 
pee pee na 25) 
ek ne Bh 
One 
ges 1 Nxt me 
J ~ aN Ty Se (26) 
d n=l 


\ 


with I, = I(k) and d the length of the segment. POA) 
memo unit pulse, nonzero only on segment n; 7 aspect bois Li 
Mector pointing in the direction of £ from one endpoint of 
Beement m to the other. The vector functions G. age 
Known as basis functions, expansion functions or dipole 
wees. ihe mode currents go to zero at the endpoints 
of the dipoles, as depicted in Figure 4. 

it is advantageous to have a symmetric impedance matrix 
Zn (to be defined later). Furthermore, the test sources 
emould be chosen to yield a well-conditioned system of 
iinear equations for the subsequent inversion process. 
For these reasons and to obtain a closed form expression 
for the integral of Equation (19), the expansion function 
for the test source Ge will be of the same form as the 
basis function for the original dipole's current distribution 


G: 


538) 








FIGURE 4. Basis functions (G,) representing the 
assumed current distributions on a 
linear filamentary antenna. 
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Winehetme DPLecewise-Ssinusoidal current distribution 
sme royeda for each segment of the antenna, then from 
Equation (22) the current on any given segment can be 
Senistaqered the sum of two sinusoidal currents, each of 
forme Vanishes at one or the other end of the segment. 
This may be visualized with the aid of Figure 4. One can 
see how the equivalent surface distribution is treated as 
a collection of (N-1) aetileraretiaye tubular dipole current 
densities on an antenna consisting of N segments. The 
Me dipole has terminal current 1 jeune eine Tune ero wayh 
the current vanishing at the far ends of the two segments 
fommane the nth junction. 

mae wsolution for the unknown complex terminal currents 
{i,t proceeds by introducing a set of (N-1) filamentary 
Gipoles each having piecewise-sinusoidal currents and 
Perm located along the axis of the original antenna. No 
Penerality is lost by assuming each test dipole has unit 
Memminal current. Reciprocity is enforced successively 
between each test dipole aaceall the tubular suriace 
Gipoles in turn. The "zero reaction condition" yields 
the following equation for the mth test dipole and the 


nth surface-current dipole (refer to Equation (15)): 


SSR ' He ds Wits ° En ds : Cai) 


@ n 


mine left-hand side of Equation (27) is nonzero only 


mien = mM. This is the result of the tangential component 
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momene clectric field on the surface of the original source 
being zero except at the delta gap. Mn which equals 

ce x f is therefore also zero unless the mth test dipole 

is co-located with the nth original dipole. Expressed 


mathematically: 


Sr -H ds. = -V fr ae § = eee 
n mi. i n pJ' mn mn n “mm mn 
n 


(28) 


The right-hand side of Equation (27) can be simplified 
by carrying out the circumferential integration of I Ea 
mean bonis is equivalent to replacing the tubular distri- 


Mmaewon with filamentary segments for calculating reaction: 


St, ee a BEXce ; E at ey een) 


n ig 


OQ, 
oP 
NH 


Equating Equations (28) and (29) one obtains 


Ws Tn o . 4 [Btn En df, (30) 


nN 


Mosumings unit amplitude for the test current and noting 


that upon co-locating m =n, 


VaSnn = 7 f Malta) Ey oT (31) 
“a 


iieminaceszration is carried out on Ts the coordinate along 


the arms of the nth surface dipole filaments. The 
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symbol Sn 1S the standard Kronecker delta, defined by 


the relations: 


= 0 m=n : (32) 
Equation (31) forms the basis of the multi-segment reaction 
feeeching technique. 


BE. MATRIX EQUATION FORMULATION 


Substitute the assumed form of Ea) (refer to Figure 


4) namely 
S iligiy vel (gee 2 ) 
— S A Giecall 
a one 7 in Pte) sin kd 
sin k(2 - £) 
A nti 
Reise erin, Son ey Sin ka (33) 


into Equation (31) and define a mutual impedance term as 


. - je ae 
: — | Be Salo Ge r= 4-1) 


mm 
ei 


a E, ae 
sin ka se 


V6 : 
. J ntl sink(Z_,, - jo) ee 


© a em ag. (34) 


c& sin kad 
n 


The dummy variable £ is the coordinate along the particular 


segment, the other variables are definedas in Equation (25). 
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ee as before, is the free space electric field intensity 
radiated by the mth test dipole with piecewise-sinusoidal 
eijerenit, speciitically for the mth test dipole composed of 
a lower segment ML and an upper segment MU, the Z component 


of the electric field is 


, 1! Saat I! 2a nw 
_ | Z, 1ML . 2M 
on Haj we Ro Mb RoMn 
\ 
= 2 LMU a “IER uy 
4 1MU OMU 
A (35) 
1MU oMU 
Where as depicted in Figure 5, Ts wr, = Is of segment 
{ = = 
ML, Ta wy I, OeecermeciueU , Ry wr, Ry of segment ML, 
Ra wy = Ry Cieasecenene MU ere... . 


Meeire 5 ao merely an extension of the single segment 
with piecewise-sinusoidal currents of Figure 2 to a dipole 
Semeosed OL two such segments. 


Equation (31) may now be written 
Vv, ore lc (36) 


where 1, is the complex terminal current at the nth junc- 
ween. the calculation are reaction between the mth 
wens Aipole and the entire equivalent surface current 
Seesivy, consisting of N-1l overlapping surface dipoles, 


feeetos Che following relation: 
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origin ML MU Zon 


FIGURE 5. The mth test dipole, composed of a lower segment 
ML and an upper segment MU 


oe 





Nicol 
i 2, eo. (37) 


Equation (37) represents one linear equation out of a 
system of (N - 1) equations that are obtained by successively 
enforcing reciprocity (or calculating reaction) between the 
(N ~ 1) surface dipoles and each of the (N ~ 1) axial test 
eeeoles in turn. The test is performed in turn between 
Pach Ot une axial dipoles and the tubular dipoles yielding 
a system of linear equations that can be described by the 


matrix relation 
7. 1) US ale ae (38) 


which is of the same form as Equation (37), except that 

ae is now a square (N ~ 1) X (N - 1) matrix and - is a 
column vector. This set of equations can be solved for 

Maes unknown junction current coefficients {1,3 Giver 2 
known voltage source is placed at the actual feed junctions 
and the voltage of the other junctions are set to zero. 

The coefficients I, of Equation (38) can be solved for 
by several techniques. Crout published a paper [7] in 
Women he discusses a method for evaluating determinants 
and solving systems of linear equations with complex 


coefficients. In the case of well conditioned matrices 


the impedance matrix can be merely inverted and multiplied 
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meee VOlLtape Matrix to yield the terminal current 


coefficients: 


=] 
[J = (2 7 [vJ (39) 


mn 


Equation (39) may be written simply as- 


~1l 
[I] = [2] [VJ (4OA) 
or equivalently: 
ere ee | Ve, (40B) 
where 
~1 
[Yj = [2] 


ime complex current coefficients I. that are determined 
from Equation (38) are used to establish the current 
distribution found from Equation (26). Once the current 
Seow ibution is known, parameters of interest such as field 
Bavverns, input Sumo cawers,  Cvc., can be calculated by 
numerically evaluating the appropriate formulas. 

A wire antenna is olouacinneye Wiel stile Wwireu ie "exc’ Led 
by a voltage source at one or more points along its length. 
Merean antenna excited at the nth junction, the applied 


Toeesace matrix of Equation (40) is 


Hy 





a 


[V4 = V | (cine) 


Oe 


meee 1S, all elements of the voltage matrix are zero except 
eiemaun, which is equal to the source voltage. The current 
distribution is established by Equation (40) which for 


the ae of Equation (41) becomes 


A 


dal 


on 


[1 
tH 
=! 
tt 
< 
<4 


(42) 


n Sil 


mn 


Note that the nth column of the admittance matrix gives 
miemecurrent distribution for a unit voltage source applied 
Geeone nth junction. iverorto mm Of une Impedance matrix 
Wiemerore gives simultaneously the current distributions 
when the antenna is excited at any arbitrary junction 
were its length. For the case where only one junction 
gs excited with a voltage Vi EMesAnpuGgeadmictance is 


peeved Girectly by the feed junction current I. as 
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x = =, (43) 


F, SOLUTION TECHNIQUES 

AS mentioned previously, integrals of the form of 
Equation (34) can be obtained in closed form for coplanar 
Seements. tin the case of non-coplanar segments, the 
Mnveprals must be calculated by numerical methods, for 
example, Simpson's rule or Gaussian quadrature. 

Upon a careful examination of Equations (34) and (35) 
[imsnould be ovbious that the integral of Equation (35) 
may be expressed in closed form with the relationship [16] 


that follows: 


Pr, exp(*jr) exp(-JjR_) 
y f ar = fexp(-jz.) ECR -z sr) 
ae. PR m m “m 
i m 
r=r 


~ exp(jz_) IIe ya 10) I ae 


cc) 


where 


Xx 
Gino) = 8 eh n(x) e+ J cos (t) ~ 1 at C45A) 
C 
x ® Ge) 
Si(x) = J aw 2, (45B) 
t | 


43 





Utilizing the relationship in Equation (44), Richmond 
and Geary published a paper [16] that developed a closed 
form expression for the mutual impedance between coplanar 
skew dipoles. However, for this antenna problem, all of 
Mmeeereactions will be performed between parallel dipoles. 
A more useful closed form expression for the problem 
Peesooing Of entirely parallel dipoles was presented by 
H. E. King [10]. The complete expression with its deriva- 
mrem may be found in Appendix A and B. Even though these 
elosed form expressions are not particularly simple, they 
feminivolve only sine and cosine integrals which have 
Perynomial algorithms; hence Zon Mies sf Ormni ss «nore 
efficient than "brute-force" integration. 

Leonard Tsai [21] formally portrayed PSRMT as a 
Peeeei cular application of Galerkin's method. Consider 


Beat the operator equation 
L(f) = ¢ (46) 


mere ide eeome Venmopecravor., is Specifically the unknown 
Serrent distribution on a filamentary dipole and the 
Sreracor L yields the fields g of a distance "a" away 
(on the original dipole). Expanding f into a linear 


@emoination of basis functions: 


Nn 


fF = dlc fF (47) 
a 


4h 





and enforcing Galerkin's method directly: 
ree Se = SE tT >... (48) 
n 


Here the a would describe the weighting functions 
corresponding to the "test dipoles" which is using the 

same expansion functions as the "original dipoles", both 
subdomain piecewise-Sinusoidal pulses. The inner product 
of the left-hand side of Equation (48) corresponds directly 


to the integral of Equation (31): 


[rt 1 eae (49) 
nN m ry 
on 


used in the reaction matching technique. Note also that 
<E5 i> = ve which is nonzero only at voltage sources. 
iimoartoliows from the "Zero reaction condition". Thus 
megemorocedures of PSRMT are those employed by a direct 


application of Galerkin's method. 
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V. DISCUSSION OF NUMERICAL TECHNIQUES AND COMPUTER PROGRAM 


Evaluation of the integrals in Equation (34) may be 
accomplished by any one of several numerical meee neneien 
Meennagues such as Romberg, Weddle, Gaussian Quadrature, 
eow, as well as by the closed form expression developed 
imeropendix A. For most integrals, no representation in 
merms Of Clementary functions is possible, and approximation 
becomes necessary. Most numerical techniques are successfully 
based upon some type of polynomial approximation: trapezoi- 
@oeerule, Simpson's rule, Romberg's method, etc. Unfortu- 
Meese chnese CLechnigues have the usual error sources such 
memeruncavion, but perhaps more interestingly they are 
harassed with the question of convergence. The question 
ieee uner as continually higher degree polynomials are 
Used or as continually smaller intervals between data 
points are used, a sequence of approximations is produced 
for which the limit of the truncation error is zero. As 
may be guessed, the smaller the interval h is made, the 
more computation time required and the greater the rounding 
Peeeecrror becomes. in other words, accuracy is limited 
for any given method by algorithm errors that ultimately 
Sescure convergence. 

itis found im practice that decreasing the interval 
between sampling points below a certain level leads to 


larger errors rather than smaller. This necessitates 
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experimentation of the various parameters available in 

Ghe technique in order to establish a trend (if one exists) 
tay gives acceptable results. Unfortunately this requires 
an assumed correct solution in order to evaluate the 
performance of the integration on a supposedly unknown 
manerral. 

Table 1 illustrates how the value for input admittance 
of a half-wave dipole varies as a function of the number 
of antenna segments and the order of the Gaussian Quadrature 
meecein evaluating the integral of Equation (34). Table 2 
is a Similar portrayal of a full-wave antenna. 

ine nerm for Table 1 and 2 was chosen to be the value of 
input admittance given by King [11]. Table 1 shows no 
convergence tendency whatsoever, while Table 2 indicates 
Cnat a lower degree polynomial may be used in evaluating 
em menue anvennia 25 adivided into more segments. The 
apparent discrepancy between Table 1 and 2 is probably 
explained by the difference in the actual current distri- 
Peet ons on the ae antennas. It is clear that the problem 
is multi-faceted. 

This poncnsence Dire bilcmers) NOUlexclusive to the 
Gaussian Quadrature technique. In fact, all of the 
numerical techniques ee troublesome parameter analysis 
gnvyolved in order to effectively use them. But a far more 
Mieervant’ and Exciting reason to avoid their use is that 
of efficiency. The closed form expression involves only 


Mm@eweovaluation of the sine and cosine integrals which have 
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ORDER OF 
GAUSSTAN QUADRATURE 





Number of Segments 


TABLE 1. 


/s 4.5 3 2.5 > 3 3. 
(Lk fs F¥ a4 3.57 ci uA 
iv £5 4 ee 45° 
16 4< Sy 657 Zi = 


Convergence analysis for the input admittance 

of a half-wave dipole. The numbers respresent 
percent variation from the King-Middleton standard 
of Y= 10.8 mmho's. (Percentages expressed to the 
nearest .5%) 
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ORDER OF 
GAUSSIAN QUADRATURE 


Number of Segments 





TABLE 2. Convergence analysis for the input admittance of a 
full-wave dipole. The numbers represent percent 
variation from the King-Middleton standard of 
y = -90 mmho's. (Percentages less than 20 expressed 
to the nearest .5%) 
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Poole aleoritvhms. Unfortunately the closed form solution 
jmeevalid in general for skewed coplanar dipoles only, 
meererore, numerical integration would be required .for 
non-coplanar elements. The same antenna problem (if 
coplanar) may be solved by the closed form technique in 

a fraction of the time that numerical methods require. 

The computer program listed in Appendix B calls upon 
Suproutine CZMNe for the evaluation of the impedance 
Weetax., itt 15 here that the method of solution may be 
Seesen. If the closed-form route is decided upon, 


(7M Cwhich is also listed) is used to evaluate a 
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VI. ANALYSIS OF RESULTS 


A. CONVERGENCE CHECK 

Figures 6 and 7 illustrate the accuracy and efficiency 
of PSRMT. Leonard Tsai [21] used Sambar iigures (the 
solution labeled "Richmond's Reaction") to show the compara- 
imme rave Of convergence of PSRMT to that of point each 
Meer criangular pulses. Note that reaction matching 
Mequiregd only six segments while for the same accuracy 
Memnco matching required five times as many segments. 
Superimposed on Tsai's comparison are those points generated 
by the closed form expressions of Equation (72) of Appendix 
‘Be 

Figures 8 and 9 illustrate the efficiency of PSRMT. 
In Figure 8, as more and more flat pulses are employed 
in Galerkin's method the more the solution corresponds to 
that of PSRMT with only six segments. Figure 9 shows how 
rapidly PSRMT relates to the moment method employing 
peecewise-sinusoidal basis functions. 
| Paeouic wl Osc xXemolities the convergence of PSRMT itself. 
The comparison is drawn between PSRMT and Hallen's equation 
employing piecewise-sinusoidal basis functions. 16 segments 
do an admirable job in relating to this half-wave dipole. 
This is sufficiently true under certain conditions (as 
will be shown) to warrant a rule of thumb: use 32 segments 


mee ewavelength. 
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Mit ommemtcea ploy Of Conductance versus susceptance 
mer thin cylindrical dipoles with radius a/a = .001588. 
fies locus spirals in as the length of the dipole increases. 
Note that for small lengths the 8 segmented solution devi- 
ates from the accepted King-Middleton value of admittance. 
Many of the plots to follow will use King-Middleton values 
for comparison, which were taken from tables published by 
R. W. P. King [11]. In the evaluation of the admittance, 
the King-Middleton procedure was again used (two iterations 
and an empirical normalization were also employed). Figure 
11 indicates the advantage of taking 32 segments per 
wavelength. 

maeyre 2 is the same type as mene ii pints tae radius 
has been increased to .00635% It is evident that for short 
merowes with large radii, the length to radius ratio for 
each segment exceeds the thin-wire approximation. That is 
iiweecnat as the length of the antenna is increased in Figure 
12, the solution becomes more convergent. 

mecures 13 and 14 portray the convergence trends of the 
Weddle numerical integration technique on a half-wave dipole. 
mmermarcated in Figure 13, the solution varies only slightiy 
with a change in the number of integration points used in 
Che integration. However, the input admittance approaches 
the King-Middleton solution most nearly when the number 
of integration points equals 29. Figure 14 shows the change 
an current distribution as the number of antenna segments 


change (with the number of integration points equal to 29). 
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Again the input admittance approaches the King-Middleton 
solution as the number of segments increase. These 

mervicular trends are vaiid only for this half-wave dipole 

of radius a/x = .001588. Other geometries must be considered 
on a one to one basis due to the behavior of Weddle's 


technique on the integral of Equation (34). 
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aeeeene RICHMOnNd's Reaction 
(6 segments) 


2) 8 eat rinacon Ss CO location 
(32 segments) 


A # & PSRMT - Closed Form 
(6 segments) 


Py ACS Pa o¥ -S 5 oF 


h/a 


FIGURE 6. Input resistance versus h/2~ for a cylindrical wire 
antenna with radius a and total length L = 2h. 
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FIGURE 7. 


h/a = 2000 





ee Richmona Ss Reaction 
(6 segments) 


29% Gor ringeon S CoO) location 
(32 segments) 


a B&B # PSRMT - Closed Form 
(6 segments 


Input reactance versus h/> for a cylindrical wire 
antenna with radius a and total length L = 2h. 
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Galerkin Method with Flat Pulses 


——_— — 6 pulses 
——— — —— |3 pulses 
20 pulses 
a——— - — — 30 pulses 
I |} (ma) 
PSRMT - Closed Form 
14 
a6 A 4 6 Segments 
Si 
(¢? Noa a eee 
_ ee 
A — \ 
i “aes ~ S x 
4 z S 
REAL SA 
$ IMAG im 





-¢ 9 


FIGURE 8. Current distribution on a half-wave dipole with 
radius a/x~ = .001588. 
(Galerkin Method with flat pulses) 
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General Moment Method 


a ae 6 basis functions 
-~ —-13 basis functions 










ZUM basis nunc tions 
- - —30 basis functions 
PSRMT - Closed Form 





& A & 6 segments 


FIGURE 9. Current distribution on a half-wave dipole with 
radius h/a = .001588. 
(General moment method — piecewise-sinusoidal pulses) 
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Hallen's Equation 








Symbo] Divisions Input Impedance 
50 86.62 + j 46.78 a 
1 | (ma) 
PSRMT ~- Closed Form 
ra) 8 BS eet A224 An 
+ 16 87.02 + j 43.45 ~ 
oe 92.25 + j 37.38 « 
lo 
$ 
b 
S 
> 
4 
of 
z 
~2 » f 6 P (.e 
z/h 


FIGURE 10. Current distribution on a half-wave dipole with 
radius a/A = .001. | . 
(Hallen's Equation - piecewise - sinusoidal - basis functions) 
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King - Middleton e 
PSRMT - Closed Form 4 8 segments 
© 32L/> Segments 


(mmho's) 






~ ip 


ee bs =] 0375 
A% tad 


; 





> 
Is” G 


(mmho 's ) 


FIGURE 11. Admittance curve for thin-cylindrical dipoles with 
radius a/x~ = .001588. 
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King - Middleton w 


PSRMT - Closed Form 4 8 segments 

+ 16 segments 

B O 24 segments 
(mmho's) *% 32 segments 





Na ¥ é | /@ se ry G 
f % . (mmho 's) 
% + 441 = 4775 
L=.6366 
a4 
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FIGURE 12. Admittance curve for thin cylindrical dipoles 
with radius a/2 = .00635. 
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BPGURE 13. 


Symbo] \ 


King-Middleton 
PSRMT NOR Ot Pus. 





Cee 4. OS cy 


(numerical integration) 5 A 9.80 - j 4.84 ~ 
Weddle's Method 17.—s«# Ono 4c 42670 > 
6 segments ) 29 = 9.53 - j 4.89 + 
4) x 9.54 -j4.71 + 

uy 
G a4 (.0 2/h 


Current distribution on a half-wave dipole with 
radius a/x = .001588 (convergence check on the 
number of integration points used in 

Weddle's technique) 
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Symbo | Y 
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I (ma) King-Middleton 8.98 - 3 4.69 «+ 
sSeaments 
/¢ PSRMT 6 A OS one ao 0 
(numerical ie ; 
integration é . Ee NS 
Weddle's method18 BH 9.14 - j 4.64 
pee ery. x 9.07 - j 4.64 
REAL =29) 
7.5 z 
S 
2.¢ 
: 7 : 
2 gf s&s $ (0 
z/h 
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FIGURE 14. Current distribution on a half-wave dipole with radius 
a/~ = .001588. 
(convergence check on the number of segments used in 
Weddle's technique) 


62 





B. ADMITTANCE CORRELATION 

Figures 15 through 18 attempt to typify the admittance 
correlation between the King-Middleton solution and various 
PSRMT solution techniques. Figure 15 and 16 depict the 
solution by Weddle's method and the Gaussian Quadrature 
method, respectively. Figures 17 and 18 can be compared 
wich Figure 11 to show the effect of changing the dipole's 
meadrus. The effect is negligible as long as the length 


memeadius ratio of each segment is maintained large. 


~ 
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FIGURE 15. Admittance correlation for thin cylindrical dipoles 
with radius a/a = .001588. (32 L/A segments) 
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FIGURE 16. Admittance correlation for thin cylindrical dipoles 
with radius a/2% = .001588. (32 L/2~ segments) 
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FIGURE 17. Admittance correlation for thin cylindrical dipoles with 
radius a/a= .003175. (32 L/~segments) 
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FIGURE 19. Admittance correlation for thin cylindrical dipoles 
with radius a/xz = .004763. (32 L/a segments) 
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fee CURRENT DISTRIBUTIONS 

macwires tout hnrough 25 present current distributions on 
Gipole antennas of various lengths and radii. On each set 
of curves is indicated the expansion set used to represent 
the unknown current. The number of segments employed in 
the PSRMT technique was based upon eeeeeme of thumb 
32 L/A 

In Figures 24 and 25, the power series [2] used as a 


femearison may be expressed as 


te) = 5 oe in 


where the Cc, s are complex constants to be determined. 
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FIGURE 19. Current distribution on a quarter-wave dipole with 
radius a/x = .001589. 
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FIGURE 20. Current distribution on a quarter-wave dipole with 
radius a/2 = .007022. 
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FIGURE 21. Current Distribution on a 3/4-wavelength 
dipole with radius a/a~a = .001588. 
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FIGURE 22. 
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PSRMT - Closed Form A 
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Current distribution on a full-wave dipole 
of radius a/z = .001588. 
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FIGURE 23. Current distribution on a 3/2 dipole of radius 
a/~ = .001588. 
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Current distribution on a half-wave dipole with 
radius a/~ = .001588. 
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FIGURE 25. Current distribution on a half-wave dipole with 
radius a/z = .005. 
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VII. CONCLUSIONS AND RECOMMENDATIONS 


A. EVALUATION OF POTENTIAL INTEGRALS 


ijcvecrals of the form 


h 70) BE) 2 
nz) = J fac = + § gicizas (5.0) 


Ry R 








where 


we) 
i 


1 \/ (z-2")° 4 ae 


Cae De + Be 


wv 
I 


are extremely difficult to evaluate for "small a" because 
Of the sharp peak in the integrand at z' = z. Although 
PSRMT has been formulated on the thin-wire assumption 
(a << 4A); for very small radius values the integrand in 
Equation (34) becomes ill-behaved due to the form of the 
Peesric i aloJke 3, (refer to Equation (35)). 

This phenomena was observed in the calculation of the 
input impedance of Figures 6 and 7. In fact, for h/~»<.3 
the value of the radius (h/a = 2000) became small enough 


momehnwart computation entirely. 
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Bee APPLICATION OF THE GENERAL PSRMT TECHNIQUE 

Although this paper dealt in particular with parallel 
dipoles only, PSRMT is useful in solving a wide variety 
of thin-wire antenna problems. With the aid of the concepts 
Peete ciprocity, reaction and surface equivalence, this 
moemmtque reduces an integral equation to a set of simul- 
Caneous linear equations. This set of equations is readily 
solved on a computer to obtain the input impedance and 
meereny AGistribution on the antenna. 

Richmond and Geary [16] presented a new form for the 
Mmmeorous near-zZone field of the linear dipole with sinusoidal 
Gyzrent distribution. This electric field expression was 
employed to derive a closed form SAG ORES for the mutual 
impedance between coplanar~skew dipoles. By employing 
Mier expression for mutual impedance, complex geometries 
and surfaces may be represented by wire grids and reaction 


pueeoreed CO yield solutions. 


C. RECOMMENDATIONS 

The theory developed in this paper is by no means 
mieertetive to the specific problem investigated. The 
data comparisons were made by employing PSRMT to equally 
Segmented, thin-wire linear, parallel dipoles. Extensions 
to this problem are obvious. First, the segments need not 
be of equal length. By merely altering the computer program, 
the specific segment lengths can be applied to the appro- 


Pioteave TOrmulas to generate a more general solution. Second, 


(7 





miemoatpoles need not be parallel. By employing Richmond 
and Geary's closed form expression, the mutual impedance 
between any two coplanar elements may be determined. 
Puuemriow Che reaction between the filamentary test dipole 
emomeme eCquivalent filamentary Cubular surface dipole must 
be examined in more detail. In the case of parallel 
elements, the reaction applied between the test dipole and 
any filamentary surface dipole is the same. But due to 
the geometry of the problem involving two non-parallel 
sepments, the mutual impedance depends upon the angular 
position of the filamentary surface dipole with respect 

to the axis. The question of how does one best react these 
eemeneos is yet to be answered. E. G. Neely will present 
Peeeeerouer's Depree Thesis to the Electrical Enginecring 
Department at the U. S. Naval Postgraduate School (in 
September 1973), in which he deals with arbitrarily bent 
monopoles above a ground plane. He is approaching this 
problem by averaging the reactions taken around the 
circumference of the tubular dipole successively. Other 


Beoreaches might be investigated. 
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APPENDIX A 


TRANSFORMATION OF EXPONENTIAL INTEGRALS 
INTO CLOSED FORM EXPRESSIONS 


iae 2 component of the electric field from a dipole 
eee a time-harmonic current distribution on each segment 
is expressed by Equation (35). For the mth test dipole, 


the dipole mode current distribution is shown in Figure 26. 


(2 








FIGURE 26. Typical mode current distribution for a 
filamentary test dipole. 
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From the expressions for the 


av the endpoints of the segments 


one obtains 


Note that 


f@emresult of substituting Equation 


eS 


t 
To wr 


t 
TOML 


t 
1) Mu 


' 
Tou 


It 


1 


i 


k(I 


k(I 


k(I 


k(I 


OML 


eML 


2MU 


OMU 


al Ose iG) -/ 


IML 


cos kd - I 


-~ I COs sd, 


cos kd - I 


IML 


1MU 


1MU 


7, 


ff 


sin 


San 


Sata 


halal 


from Figure 26, it is apparent that 


IML 


oML 


OMU 


t 
To wt 


t 
Tomy 


=e 


oMU 


eae 


IMU 


oy Soa ik 


KGocmedmy Sin ikd 
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kd 


kd 


kd 


kd 


(32) into Houation 


G€erivatives of the current 


(iene Bouation (24) ),; 


(51A) 


(5a 3) 


Gre) 


Gry) 


(52A) 


(52B) 


(52C) 


al) 


(53A) 


(53B) 





mecetine the distance designations of Figure 5 according 


ier raigure ¢/. 
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Equivalent Surface Dipole 
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FIGURE 27. Geometry of two parallel dipoles in echelon. 
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Vauhieiie —feomeury indicaved in Figure 2/ and Equation 


(53) the electric field of Equation (35) may be rewritten 


as 
-jkR -jkR 
= Z ke . pee : 
m Tajwe Ry Sin ka Ro Sane kd 
-JkRo 
= \ 
Ele cos. kd R) sin kd (54) 
peecauivalently: 
i r -JkR, . -JkR, 
EF = Hitt ——— ) o — OO he aa Oe aren Aa 
m nw ¢ R, sin kd R, Sin kd 
-jkR 
= 0 
+ 2k cos ka poems namical (55), 
in this problem 
k = 27 /AaA with OMe ee eS al 
w= 2n(3(10°)) 
2S 1/aeae (or 
therefore 
k/utwe = 30. (56) 
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Eeametion (55) in light of Equation (56) is 


7 as ‘ ao ese 2 26 
ee elke eR. CS COS Ke 
al 2 0 
rat) 
where 
Ro = ae 4. z° (58A) 
_ 2 2 

R, == a + (Ly ~ Zz) (58B) 


R, = ae (2, + 2)° (58C) 


Piproying the trigonometric identity 
See = ce Sins COS <<  , 


Equation (57) may be written as 


-jkR, -jkR., ~JkRy 


fe 20 2 -~j e hy aie A comeosnkd Ze 
m sin kd Ry Ro Ro 


yD) 


Equation (59) is exactly the form that H. E. King expressed 
ieee lectric field in; except for the extra’ sin ka’ term ng) 
the denominator. This term is not present in King's 


expression for the electric field because he was assuming 
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beateewne anvenna current distribution is sinusoidal. 
ioeene application of PSRMT the current distributions are 
assumed to be piecewise-sinusoidal and therefore the 


expression for the mutual impedance is: 


L. | 
a = 30 __, | eaten (2 cae) 
le (sin ka) g i 


n-l 
ie 
f ; sin Kner 
n 
-JkR, -JkR, -JkR, 
-j e ff -~je A ZL Cos kde ar 
R Ig: R 
1 2 0 
(60) 


memaeion (60) follows from the substitution of Equation (59) 
ice Equation (34) and assuming that the length of all 
segments are equal to d. 

imeorder to facilitate transforming Equation (60) into 
mmetosed form expression, Zon will be expressed in terms 
oeome feometry of Figure 2) Omens oldahaourrent 


Sewer lOutTLONS, 


-— 3 f 2k, th 
Zn = naa J sin k(zg —-h) + Sn) K(2tytn~2)} 


£,+h 
~jkR -jkKR -jkR 
-je - -j e a 2] cos k Gl 0 aver 
Ry Ry Ro 
(61) 
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Obviously Equation (61) can be multiplied out to yield 


pemeee xpressions, one of which follows: 


ek,th -~jkR 
a = +30 Salial k (2k, +h - Zz) pp lik Aitreeteen az. 


mipeeying the trigonometric relations: 


e = Wises CS eo liga (63A) 
Sins = -sin (-) (63B) 
cos“* = cos (-«) , (63C) 


Equation (62) becomes 


| Sone 1 COS KR, Sain 
ih 
Leth 1 il 


a = sin k(2£L,th-z )dz 


(64) 


The real part of complex Equation (61) gives the mutual 
mectovance and the imaginary part gives the mutual reactance. 
The real part of Equation (64) gives only one-sixth of the 
Pele expression for the mutual resistance, but it will now 
Bemeranstormed into a closed form expression in order to 
illustrate the procedure of obtaining all twelve Terms , 


both resistive and reactive. 
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| ool ae Pe eer 


R 
L,+h it 


A = 30 


REAL Sime 22a he 2) az 


co) 


Substitute the expression for R, given in Equation (58) 


al 
momooptain 


| Sigs sin k ae + ( a z)° 
Appar, a a 5 ; sin k(2£,th-z) az. 
2 a + ( Sy D)) 
(66) 
By a change in the variable of x = tL - Z, Equation 


(66) becomes 


L,-2L,—-h aes 
A = -30 | sin kVA + x 
REAL a a 


sin k(2k,th-£, +x) eer 


(67) 


Employing the trigonometric identity for the product of 


angles one obtains 
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£,-2£,-h 
Tel) f [-cos k(h-£, +22, ) cos k( ActN 2 
£,-£,-h 


2 


- sin k(h - £, + 2k.) Sine ae thse -~ xX) 


+ cos k(h - fy + 2k.) cos k( ae bh xX) 


-~ sin k(h - Ly + 2k.) Salat Kat tx + x)] > ee cesas 
ay) 2 ® 
a ex 
(68) 
feoecond change of variable of 
wl ee let atx ae) (69A) 
y =k aoe + x) (69B) 
gives 
3 
x cos u 
Appar, = sf [+cos k(h - £, i 2k.) ras au 
1 
: , Sailiey wl 
+sin k(h - Ly SE el.) aac: du 
te 
+ 15 f [+eos k(h - £, + 2k.) ———. Gv 
5) 
Dey _ Sane Vv: 
sin k(h £, + 2k.) cae av eee 
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where 


7 wie 2 
Oo ees et Ch — £, + £,) + (h - £, + £,)) 
nee k( Vac + (h - £, + 22 \e + (h =f) ee) ) 
3 l Z al 2 
: V2 2 
a = mera +: = £, + £,) oh é 4 £,)) 
_ = ieee. 2 
V2 = el ee 2 one a = 2L.) ~ (h - £, + 2k,)) 


The integrals of Equation (70) are recognized as sine 


and cosine integrals and may be written as 


Memaz, = 15 | cos k(L, = ans = he fin, = Ci(v,) + Ci(uz) + Ci(v,)] 


ar lal! KCL, - ek. - h) (Siu, ) ~ Si(v,) - Si(u,) a Si(v3)] 


a 


femal that Equation (61) through Equation (70) was written 
mest nusoldal current distributions. To make them appli- 
cable to piecewise-sinusoidal current distributions, 
Homerton (61) must be multiplied by the factor l1/ (sin ce 

Equation (61) expanded in closed form as typified by 
Beton (71) is written in its entirety in Equation (72) 
and (73) in Appendix B. Equation (73) is an expression 


for the mutual impedance between two thin parallel center-fed 
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paeennas With half-lengths HL spaced a distance DL apart 
and staggered by the value SL. This expression is a 
emolsttcation of H. E. King's general formula for mutual 


impedance [10]. 
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APPENDIX B 


FORTRAN PROGRAM LISTING 


Pigeure 28 is the main computer program listing (in 
memeran 1V) for the solution of the linear dipole antenna's 
Seevine point admittance and current distribution. 


dime symbols for the input data are defined as: 


EL length of the dipole in wavelengths 
A radius of the dipole 


N number of antenna segments. 


imis program has built in “"do-loop" that generates 
mieevoltage matrix of Equation (41). The voltage source 
Semone VOlt is automatically placed at the center junction 
Oi the dipole. 

Calculation of the complex ae matrix is done in the 
Suietoeucine CZMN2. This calculation may be performed by 
@m@emo: Several techniques. Gaussian Quadrature, Komberg 
and Weddle's methods can be employed in order to effect 
numerical integration of Equation (34) for es nf wads 
memoone, additional subroutines are required to calculate 
meemeorectric field explicitly, establish the limits of 


Mame ration, perform the necessary integration technique, 


a 





me .. cae may be most efficiently calculated by programming 
directly the closed form expression of Equation (73). 

Once re was determined, CMIN1, a subroutine developed 
Myers Gc. Neely, was used to calculate the complex inverse 
One A ine complex current coefficients result from a 
direct application of Equation (40). Since a single one 
Weeeesource was placed at the center junction, the driving 
point admittance is known immediately from Equation (43). 

The program's output is set up to provide a.) an echo 
eheck on the input parameters and important variables, 

b.) list the real and imaginary parts of the complex 
[miecion currents as a function of z/h, c.) list the 
magnitude and phase of junction currents, am@eed.) display 


the input admittance and impedance. This output is 


presented for both the dipole and the monopole. 
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